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ON  SUMDER  B(»1ES  OF  MINIHUM  DRAG 
IN  NEWTONIAN  FLOW 
by 

DAVID  Q.  HDLL^*^ 


SOTftART 

In  recent  papers  by  Mlele  (Refs*  1  through  3)«  the  problea  of  sdnl- 
■Izlng  the  pressure  drag  of  slender  bodies  in  Newtonian  flow  ws  considered 
in  general  and,  then,  solved  for  those  particular  oases  in  which,  of  the 
four  geoowtrie  properties  being  considered  (thickness,  length,  enclosed  area, 
and  noaent  of  inertia  of  the  contour  for  two-diaensional  shapes  and  diaaeter, 
length,  wetted  area,  and  volume  for  axisynmetrlc  shapes),  two  are  prescribed 
and  the  reaainlng  two  are  free*  In  this  paper,  the  analysis  of  Refe*  1 
through  3  is  extended  to  the  class  of  problems  in  which  three  quantities  are 
prescribed  and  the  remaining  is  free*  After  the  variational  problea  is  re- 
foranilated  in  order  to  account  for  the  fact  that  the  pressure  ooefficieat  must 
be  aonnegatlve  everywhere,  special  attention  is  devoted  to  those  particular 
cases  in  which  two  of  the  three  prescribed  quantities  are  the  thickness  and 
the  length*  In  each  case,  a  one-parameter  family  of  extreaal  solutions  is 
obtained,  the  parameter  being  related  to  the  three  prescribed  quantities* 
Furtheraore,  each  family  of  extreaal  solutions  contains  three  classes  of 
body  shapesi  (I)  an  infinitely  thin  plate  or  a  spike  followed  by  a  regular 
shape,  (II)  a  regular  shape  only,  and  (III)  a  regular  shape  followed  by  a 
constant  thickness  contour  or  a  cylinder*  In  all  of  the  eases  considered, 

^*^Staff  Aesoeiate,  Astrodynaaies  and  Flight  Mechanics  Group* 
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analytleal  •xprtaslona  art  obtainad  for  tha  gaoaatry  of  tha  optiaua  ahapas 
and  tha  aaaooiatad  drag  eoaffieianta* 
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1*  IWTROPqCTlOlt 

Zo  tlw  preTloua  raporta  by  Miela  (Rafa*  1  and  2),  tha  problaa  of  aial* 
■ising  tha  praaaura  drag  of  alaadar  bodiaa  in  Nawtooian  flow  waa  fonnilatad 
la  gaaaral  and  aolwad  for  thoaa  paurtioular  eaaaa  in  tdxleh*  of  tha  four  gao« 
■atrio  propartiaa  balng  eonaidarad  (thloknaaa  d«  langth  i|  anoloaad  araa  K, 
and  aoaant  of  inartla  of  tha  contour  N  for  two-dimanalonal  ahapaa  and 
thlcknaaa  d,  langth  i,  %fattad  araa  S*  and  voluaa  V  for  axlayanatric  idiapaa), 
two  ara  praaeribad  and  tha  raoadnlng  two  ara  ftaa.  In  thla  raport,  tha 
analyala  of  Rafa*  1  and  2  la  axtandad  In  ordar  to  cover  tha  caaa  in  «diieh 
three  of  tha  gaoaatrio  propartiaa  ara  apacifiad  and  only  one  ia  Area* 

However,  in  ordar  to  do  ao,  it  la  nacaaaary  to  raforanilata  tha  variational 
problaa  and  include  tha  condition  that  tha  praaaura  eoafflciant  auat  be 
nonnagativa  avarywhara.  In  tha  intaraat  of  brevity,  tha  approach  of  Raf*  3 
ia  aajdloyad,  that  is,  tha  two>diaansional  problaa  and  tha  axisyaBatrie 
problaa  ara  eonaidarad  siaultaneoualy*  In  particular,  attention  ia  devoted 
to  thoaa  oaaaa  in  which  two  of  the  three  praaeribad  quantitiea  ara  tha  thiek- 
naaa  and  tha  length. 
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2.  MnnmiM  drag  wotiim 

Aoeordlag  to  Rofa.  1  through  3%  tho  Noiftonian  preseurt  ooofficlont  for 
a  alandor  body  in  an  Invlaold  hyparaonlo  flow  at  aero  angle  of  attack  la 
given  by 


C  -  2ir  Cl) 

P 

where  x  denotes  a  coordinate  in  the  flow  direction ^  y  a  coordinate  noraal 
to  the  flow  direction «  and  y  the  derivative  dy/dx*  With  reference  to  the 
portion  of  the  body  between  stations  0  and  x,  the  drag  per  unit  spant  the 
enclosed  area,  and  the  aoment  of  inertia  of  the  contour  of  a  two-dlaeneional 
shape  and  the  drag,  the  wetted  area,  and  the  voluae  of  an  axisyaaetrio 


shape  can  be  written  as 

Twoi«diaeneional  case 

D(*)  -  4q  dx 

A(x)  -  2  y  dx 
P*  2 

M(x)  -  2  y  dx 
Jo 

After  the  definitions 

Two-dineneional  case 

or  -  D(x)Aq 

9  -  A(x)/2 


Axisyaaetrio  case 
D(x)  ■  41*1  yy^  dx 
S(x)  ■  2n  (2) 

V(x)  ■  n  J*  y^  dx 

Axisyaaetric  case 
a  ■  D(x)/4Tai 

9  ■  S(x)/2rt  (3) 


Y  •  M(x)/2 


Y  ■  V(x)/n 
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•re  latreduoed  and  after  both  aides  of  Eqs.  (2)  are  differentiated  with 
respect  to  the  independent  variable)  the  following  differential  oonstraiats 
are  obtained  (Ref*  3)i 


•  n«3  _ 

«  -  y  r  ■  0 

0  >  y  a  0  (^) 

•  2  _ 

V  -  y  *  0 

where  a  «  0  in  the  two-dinensional  ease  and  n  ■  1  in  the  axiayometrio  oass* 
Since  the  condition  that  the  pressure  coefficient  mist  be  nonnegative  can 
be  expressed  in  the  fom 


y  -  p^  »  0  (5) 

where  p  denotes  a  real  variable)  the  systea  of  differential  equations  (4) 

and  (3)  has  one  independent  variable  (x))  five  dependent  variables 

(cTi  0)  Y*  7f  P)(  degree  of  freedoa*  Consequently)  if  it  is  assuned 

that 


Xj^  ■  y^  a  a  0J^  »  ■  0  (6 

and  that  soas)  but  not  all)  of  the  reaaining  state  variables  are  specified 
at  the  final  point*  the  mlniaua  drag  problea  is  stated  as  followst  In  the 
class  of  functions  oKx).  B(x).  v(x).  y(x))  p(x)  which  are  consistent  with 
the  differential  constraints  (4)  and  (5)  and  the  initial  conditions  (6)« 
find  that  special  set  which  ainiaizes  the  difference  IB  •  where 


a  ■  or* 
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3.  HECESSART  COMDITIONS 

This  problem  ia  one  of  the  Mayer  type  with  separated  end-oondltlonst 
80  that,  after  the  Lagrange  multlpllera  through  Xj^  are  introduced,  the 
extremal  arc  Is  governed  by  the  following  Euler-Iagrange  equations  (Refs*  8 
and  9): 

Xi-O 
X2  »  0 

Xj  -  0  (7) 

^  (-  3X^y”y^  +  X4)  -  -  Xiny”“V  -  X2  -  ZX^y 

0  -  -  2X4P 

Integration  of  the  first  three  eqtiatlons  leads  to  the  results 
Xj^  ■  Cl  ,  ^2  “  S  “  ®3 

e 

where  C^,  C2,  are  constants*  Furthermore,  since  the  independent  variable 
does  not  appeeu:  explicitly  in  the  constraining  equations,  the  additional 
first  integral 

-  2Cj^y“^  +  C2y  +  Cjy^  ■  C  (9) 

is  valid,  where  C  is  a  constant* 
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Ab  th*  fifth  Eulsr  •quation  indicates,  the  extremal  arc  is  generalljr 
discontinuous  and  is  coaqwsed  of  the  suboros 

•  0  and/or  p  ■  0  (10) 

Along  the  former  subarcs ,  called  regular  shapes,  ths  pressure  coefficient  is 
alMsys  positive  or,  at  most,  zero  at  a  single  point*  Along  the  latter  sub- 
arcs,  the  slope  of  the  body  is  always  zero;  hence,  they  are  called  zero- 
slope  shapes* 

When  a  discontinuity  occurs,  the  Erdmann-Welerstrass  comer  conditions 
must  be  applied.  They  require  that  the  constants  C^,  C^t  C  have  ths 
sane  value  for  all  of  the  subarcs  included  in  the  extremal  arc  and  that 

A(-  3X^A*^  ♦  \^^)  •  ■  0  (11) 

Consequently,  Sqs*  (11)  imply  that  ■  0  at  a  corner  point*  Furthermore, 
a  discontinuity  in  the  slope  is  not  possible  in  the  two-diaensional  case 
(n  ■  0)  but  is  possible  in  the  oxisynmetrlo  case  (n  >  1)  if  y  ■  0* 

The  end-conditions  are  partly  of  a  fixed  end-point  type  and  pertly  of 
the  natwol  type*  The  latter  must  be  derived  from  the  tronsversality  con¬ 
dition 

Cdx  +  (Oj^  >  1)  da  -I-  Cgdp  +  Cjdy  +  (-  ♦  X|^)  dyj^  ■  0  (12) 

which  is  to  be  satisfied  for  all  systems  of  differentials  consistent  with 
the  prescribed  end-conditions  and  implies  that  •  -  1* 
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Tb*  applieatioB  of  tho  Logondre-Clobach  condition  indicates  that  the 
drag  is  a  sinimB  if  the  following  inequalities  are  satisfied t 

j  >  0  t  along  the  regular  shape 

(13) 

<  0  t  along  the  sero-slope  shape 

Froa  the  previous  dlacussion,  it  appears  that  the  multiplier  plays 
an  laq[)ortaat  role  in  determining  the  composition  of  the  extremal  are  and« 
hence «  is  called  the  switching  function*  Its  properties  are  summarized  as 
follows! 

X^  ■  0  t  along  a  regular  shape 

X|^  ^  0  «  along  a  zero»slope  shape  (l4) 

X|^  •  0  t  at  a  comer  point 
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qPnPAL  flOLBPIOW 

Aft«r  introduoing  th*  noBdimtimional  vorlAblca 


«-7 


(15) 


«Bd  tb*  dcfiiiltlo&a 


-  C(d/2)’  “  ,  *2  -  C2(d/2)^  “  "  ,  Kj  -  Cj(d/2)^  "  “  (16) 


OM  can  rtMrit*  th«  firat  iatagral  (9)  in  tha  fom 


(17) 


la  which  t  ■  d/i  daaotaa  tha  thioknaaa  ratio*  Dua  to  tha  fact  that  tha 
praaaura  ooaffioiant  nuat  ba  noanagatiaa  ararywharaf  tha  applloatioa  of 
thla  firat  iatagral  at  both  aad  pointa  of  tha  axtraaal  arc  laada  to  tha 
baaic  iaaqualitiaa 


K^aO,  Kj-Kg-KjaO 


(18) 


rurtharaora*  tha  gaaaral  aquatioaa  for  tha  regular  ahapa  aad  tha  aaro- 
al<^  ahapa  can  ba  written  aa 


-  K-  l^ 


dTt  ♦  conat 


Kgll-  Kj 


(19) 


■  Conat 


raapaotiaaljr* 
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If  d«Bot«a  th*  ordinate  of  a  oomer  point  between  a  nero-alope 
shape  and  a  regular  shapet  the  transition  eqmtion 

'2’!.  *  'S’€  ■  *1  <“> 

holds*  Furthemore,  after  the  Bondimensional  switching  function  is  defined 
as 

(a) 

the  integration  of  the  fourth  Euler-Ingrange  equation  along  a  sero-slope 
shape  yields  the  relationship 

O  •  K(§g  -  f)  ,  K  -  Kg  +  (22) 

This  equation  supplies  the  proper  sequence  of  subaros  as  foUowsi  (a)  If 
K  >  0,  a  regular  shape  precedes  the  comer  point*  and  a  zero^slope  shape 
followsi  furthemore,  owing  to  the  monotonic  nature  of  the  switching  function* 
no  subaro  nay  follow  the  zero-slope  shape*  (b)  If  K  <  0*  a  zero-slope 
shape  precedes  the  corner  point*  and  a  regular  shape  follows|  also,  no  sub- 
arc  nay  precede  the  zero-slope  shape. 

In  view  of  these  stateaents  and  the  transition  equation*  the  nairinun 
possible  nuaber  of  subaros  coaposing  the  extreaal  are  is  three  (one  regular 
shape  and  two  zero-slope  shapes)  with  one  of  the  two  possible  corner  points 
lying  along  1)  ■  0  and  the  other*  along  T|  ■  1*  If  the  coordinates  of  these 
two  comer  points  are  denoted  by  ^*  0  and  ^*  1*  the  general  equations  for 
the  shape  of  the  optiaua  body  can  be  expressed  as 


u 


.  t| 


dll 


5-?o 


r*/zzzzz 

Jq  Jk^  -  KgD  -  IC,Tp 

r 

Jo  7Ki-K2T|-lC3r 


(23) 


?  <1  .  D 


For  oxtroaal  ores  involving  onlj  ont  comer  point  or  no  corner  points  at 
all,  the  dliaps  of  the  optiaun  body  can  be  formally  obtained  from  Bqs*  (23) 
by  setting  ^  ■  0  an4/or  ^  •  1,  It  must  be  understood,  however,  that  the 
comer  conditions  need  not  be  satisfied  at  these  points* 

Once  the  shape  is  known,  the  following  nondimensional  integrals  can  be 
evaluated! 

I„  .  Tft?d?  .  W5  ,  I„  .  I,  .  ifds  (2*) 


SO  that,  after  the  corresponding  dimensional  quantities  are  written  as 


Two-dimensional  case 

Axisymmetric  ease 

2r  " 

«nd'» 

A-dfl^ 

S  ■  tidje  Ig 

M 

(25) 
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It  is  peasibls  to  oxproao  ovary  uakaowii  quantity  in  torna  of  tha  kaowi 
quaatltiaa  and  tha  nondlnanalonal  intagrala  (24)*  finally,  If  tha  rafar* 
anca  araa  for  tha  drag  eoaffielant  la  ehoaan  to  ba  tha  frontal  araa  avalnatad 
at  X  ■  it  tha  following  ralationahip  can  ba  aatabll^ad  batwaan  tha  drag 
eoaffielant  and  tha  thlcknaaa  ratio t 


(26) 


a 
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3*  OPTIMai  l!WM)lMENSIOmL  BODIES 

In  th*  prtvloua  Mctlona,  th*  alalmai  drag  problM  was  aolTad  in  gtaaml 
for  arbitrary  boundary  conditions*  Bars*  two  particular  casts  associatsd 
with  ths  two>diMasioaal  problwa  (a  •  0)  are  analysed  in  detail^  and  the 
results  are  sunoairlsed  in  Figs*  1  through  7* 

5*1*  Given  Thickness*  Lenitth*  and  Enclosed  Area 

If  the  thickness,  ths  length,  and  the  enclosed  area  are  prescribed 
lAUe  the  anaeat  of  inertia  of  the  contour  is  free,  the  transvarsality  con¬ 
dition  leads  to  ■  0  which  inplies  that  ■  0,  so  that  the  basic  in¬ 
equalities  (l8),  the  transition  equation  (20),  and  the  witching  function 
(22)  reduce  to 


K^aO,  K^-KghO 


1]^  -  KjA2  #  a  -  -  5) 


(27) 


These  results  indicate  that  the  comer  point  0  exists  tdien  1^-0,  which 
inplies  that  <  0*  On  the  other  band,  the  comer  point  1  exists  whsn 
■  I^,  which  inplies  that  >  0*  Consequently,  these  two  comer  points 
cannot  coexist*  lUrthernore,  for  >  0  and  >  0,  no  comer  points 

exist* 

In  conclusion,  three  classes  of  bodies  enter  into  the  coe-paraneter 
fanily  of  extrenal  solutions  (Fig*  l)t  (I)  an  infinitely  thin  flat  plate 
followed  by  a  regular  shape,  (II)  a  regular  shape  only,  and  (111)  a  regular 
shape  followed  by  a  constant  thickness  contour*  The  paraneter  which  governs 
these  solutions,  called  the  shape  paraneter.  can  be  defined  as 
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Ip  ■  ^  (28) 

Bodies  of  Claw  1«  This  elsss  of  bodies  •  which  is  ooaposed  of  in- 
finitelgr  thin  flat  plates  followed  by  regular  shapes,  is  characterized  by 
the  conditions 


Kj-O,  K2<0,  0s?^<1.  5^-1  (29) 

which  are  valid  for 


0  <  «p  s  I  (30) 

Consequently,  use  of  Eq*  (23)  yields  the  following  expression  for  the  shape 
of  the  optisms  body  (Fig.  2): 

0  <  5  <  .  1)  -  0 

/f  - 

where,  fron  Gq.  (28),  it  is  seen  that  the  abscieea  of  the  transition  point 
is  related  to  the  shape  parameter  by  (Fig.  3) 


5,  •  1  -  ?  02) 

Finally,  froa  Eq.  (26),  the  sinisus  drag  coefficient  can  be  written  in- the 
toxm  (Fig.  4) 
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T  250  (jT 


(53) 


Bodlf  ot  Cl*—  n*  ni»  mxt  olkM  of  bodioo  ia  oe^toaod  of  rogalar 
ahapoa  only  and  ia  oharaotarlaad  by  th#  oonditloaa 


K^>0*  K^-K2>0,  ?^-0, 


(34) 


whieh  hold  for 


I  a  (p  a 


(35) 


nuMf  th#  ahapo  of  tha  optlnaa  body  la  glvon  by 


?/5 


0  a  ?  al  ,  5  , 

1  .  (1  - 


(36) 


tdiora 


t  ■  K^Aj^  t  -  «  a  ♦  a  1 


(37) 


ia  a  ooaataat  whieh  ia  related  to  the  ahape  paraaeter  through  the  expreaaloB 


. .  >  ?  -  g  ♦  it? 

’  ^  1  -  (1 .  ,)^ 


(36) 


The  drag  ooefflolent  can  then  be  written  aa  foUowat 
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(39) 


Bodlea  of  Claae  Hit  Th«  last  class  of  bodies  is  ooapossd  of  regular 
shapes  followed  by  constant  thickness  contours  and  satisfies  the  oondltlons 

ICj>0,  Sq-O.  (40) 

which  are  valid  for 


<  cp  <  1 


(41) 


The  equation  for  the  shape  of  the  optima  body  is  then  given  by 


<fslt  Tl“l 


/2 


(42) 


idiere  the  abscissa  of  the  transition  point  and  the  shape  paraneter  are  re¬ 
lated  as  follows I 


■  I  (1  -  «p) 


(43) 


Ths  associated  drag  coefficient  can  be  expressed  in  the  form 


•r  250(1  -  cp)‘ 


JL 


(44) 
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Th*  above  roBUlta  eaa  bo  o^pleyod  in  order  to  understand  the  effect 
of  changing  any  one  of  the  given  quantities  on  the  shape  of  the  optlmn 
body  and  Its  drag  while  the  other  two  are  kept  constant*  As  an  ewqilet 
If  the  thlckneas  and  the  length  are  constant  and  the  enclosed  area  Is 
varied,  then  the  pressure  drag  (Fig.  b)  is  an  absolute  ninlnwi  when  f  ■  1/2. 
The  corresponding  optinsa  shape  is  a  wedge,  which  is  precisely  the  result 
obtained  by  Miele  in  Ref.  1  for  problems  where  the  thickness  and  the  length 
are  given  while  the  enclosed  area  and  the  aoaent  of  Inertia  of  the  contour 
are  free* 


9.2.  Given  Thickness.  Length,  and  Moment  of  Inertia  of  the  Contour 

If  the  thickness,  the  length,  and  the  moment  of  Inertia  of  the  contour 
are  prescribed  while  the  enclosed  area  is  free,  the  transversality  oondition 
leads  to  C2  ■  0  iMoh  means  that  K2  ■  0.  Thus,  the  basic  inequalities,  the 
transition  equation,  and  the  switohing  function  become 


,  o  -  -  ?) 


(‘»5) 


and  imply  the  foUowingt  (a)  The  corner  point  0  occurs  when  ■  Of 
(b)  The  comer  point  1  occurs  when  ■  K^t  and  (c)  These  corner  points 
cannot  coexist* 

Again,  one  obtains  a  one-parameter  family  of  extremal  solutions  which 
is  composed  of  the  three  classes  of  shapes  mentioned  in  the  previous  seotion 
(Fig.  1).  The  shape  parameter,  which  governs  these  solutions,  is  defined  as 


Bodiw  of  Claea  I«  Th«  conditions  to  bo  oatlofied  for  this  class  of 


bodies  are 


K^-O,  1C3<0,  0<C^<1, 


(•*7) 


and  are  valid  for 


0  <  «p  <  i  (48) 

Consequently,  Eq*  (23)  yields  the  following  expression  for  the  optimn 
shape  (Fig*  5)t 


0  «  ?  s  . 

«  1  . 


-  0 


(49) 


where  the  absolssa  of  the  transition  point  is  related  to  the  shape  paraneter 


by  (Fig.  6) 


-  1  -  7® 


(50) 


Furthernore,  the  mlnianin  drag  coefficient  is  given  by  (Fig.  7) 


s 


“2 


(51) 
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Bodif  of  Cl*—  II.  This  class  of  bodies  is  eharsetsrissd  bf  ths  oob- 

ditiOBS 


lt^>0,  Kj-ICj>0.  ?^-0.  ?^-l 


«diieh  arc  valid  for 


1 


7 


a  Ip  < 


(52) 


(53) 


IhttSt  ths  sxprssaioa  for  ths  optiaua  shape  is  as  follows! 


idlers 


♦ 


« 


K^<0  Osts* 

Kj  »  0  0  «  ♦  <1 


sod  where 


A(1|,  «)  -  i 


J*  tl 


♦  F(9,  k)  -  K(t, 


k) 


(5*») 


(55) 


(56) 


la  these  aad  the  subsequeot  foraulas  of  this  seetioa,  the  upper  sigas  are 
to  be  si^oyed  idiea  s  0  which  aeans  that  s  ip  a  1/3*  and  the  lower 
sigast  whea  h  0  which  aeaas  that  1/3  a  ip  <  3/^*  Ttartheraore*  the  syibbolB 
T  aad  I  deaote  the  iaooaplete  elliptic  iategrals  of  the  first  aad  second 
kiad|  respectively,  whose  argowat  9  aad  paraaeter  k  are  defined  as 
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e  -  arc  eo8  ^  ^  (?7) 

‘■7^ 


Next,  tha  conataat  #  la  related  to  the  shape  paraneter  throu^  the  expression 


(58) 


tdiere 


B(T|,  1K1  *  tlf 


(59) 


Finally,  the  alnlaun  drag  coefficient  la  given  by 


7 


(60) 


Bodies  of  Claaa  1II»  The  bodies  of  this  class  oust  satisfy  the  follow¬ 
ing  conditions: 


Kj^>0,  K3-K3,,  0<?^<1 


(61) 


which  are  valid  for 


St  cp  <  1 


(62) 


I 


r; 


n 


w 
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CoiiMqu0iitly«  th*  shtp*  of  th«  optinai  bodjr  can  ba  axpraasad  as 


(63) 


?,<?<!,  1)  -  1 


whtra  tha  abaolaaa  of  tha  transition  point  is  givan  by 


^  (1  -  <p) 


(64) 


Tha  drag  coaffioient  for  the  optiaua  body  cam  than  be  witten  as 


(65) 


i 


I 
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6.  OFTIWM  AXlSTfBTBIC  BODIBS 

m  this  ssetioBi  two  particular  oases  associated  with  the  three- 
diaeaslonal  probles  (n  -  1)  are  oonsidered,  and  the  results  are  susnarlsed 
in  Figs*  8  through  13* 

Oiren  Thickneaa.  Length,  and  Wetted  Area 

If  the  thickness t  the  length*  and  the  wetted  area  are  given  while  the 
volune  is  free*  the  transversallty  condition  leads  to  ■  0  and  isplies 
that  ■  0*  ao  that  the  basic  inequalities*  the  transition  equation*  and 
the  switching  function  reduce  to  expressions  (27)*  Consequently*  the 
comer  point  ^^*  0  exists  for  ■  Otand  the  corner  point  Cj^*  1  exists  tdwn 
■■  Incidentally*  these  two  corner  points  cannot  coexist*  Further- 
■ore*  for  >  0  and  >  K2  >  0*  no  corner  points  exist  (Fig*  1)* 

It  is  evident  froa  the  previous  discussion  that  three  classes  of  bodies 
enter  into  the  ons-paraaeter  faaily  of  extreaal  solutions t  (I)  a  qpike 
followed  by  a  regular  shape*  (II)  a  regular  shape  only*  and  (III)  a  regular 
shape  followed  by  a  cylinder*  The  shape  paraaeter  governing  the  solutions 
is  defined  as 


Bodies  of  Class  1*  For  this  class  of  bodies*  \idiich  is  coaposed  of 
spikes  followed  by  regular  shapes*  conditions  (29)  aust  be  satisfied  and  aure 
valid  for 


0  <  tp  s 


1 


(67) 
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ConMqiMBtlyt  um  of  Eq*  (23)  ylolda  tho  foUowliig  roault  for  tho  optlMai 
■bopo  (llg*  8)i 


(€8) 


Tho  absolaoo  of  tho  traBoltioB  point  can  bo  obtolnod  froa  Eq.  (66)  OBd  !• 
«lvoB  by  (Fig.  9) 


^  .  1  -  2v 


(69) 


flBollyt  tho  BinlwiB  drag  cooffioiont  can  bo  exprtaoodas  (Fig.  10) 

^  -  -ij  (70) 

T  8«p 

Bodioo  of  Claoo  II.  Thlo  olaoo  of  bodloa*  whioh  io  ooqtoood  of  rog«> 
lar  ahapos  oaly,  la  obaraotoxlsod  by  tho  aot  of  oonditioBa  (3^)  idiioh  aro 
valid  for 


1  2 

Tho  oiqiroaaioii  for  tho  ohapo  of  tho  optiana  body  ia  givoB  by 


(71) 


0  <  ?  d  1  , 


(72) 
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whert 


♦  ■  *2^^  •  -  •  <  ♦  <  1 


and  whar* 


A<tl.  «) 


+ 


i  are  tan 

S 


■■  .AM _ 


(73) 


The  constant  f  is  related  to  the  shape  paranster  throuj^  the  expression 


«P 


tdiere 


B(T.  t)  •  (1  - 


The  assoolated  drag  coefficient  can  be  expressed  as 


(74) 


(75) 


(76) 


Bodies  of  Class  III«  The  set  of  conditions  to  be  satisfied  by  this 


class  of  bodies,  which  is  eonposed  of  regular  shapes  followed  by  cylinders, 
is  represented  by  Eqs*  (40)  and  hold  for 
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I  <  «p  <  1  (77) 

Tlw  equation  for  the  optiaua  shape  can  then  be  written  ae  follows: 

(78) 

where  ths  absolssa  of  the  transition  point  is  given  by 

-  3(1  -  «p)  (79) 

rinalljt  ths  drag  coefficient  can  be  expressed  as 


6,2,  given  Ihicknssst  Length,  and  Volume 

Xf  ths  thickness,  the  length, and  the  voluaw  are  given  while  ths  wetted 
area  is  fk^s,  the  transvsrsadi^  condition  leads  to  C2  ■  0  and  ijg>lies  that 
K2  ■  0,  so  that  ths  basic  inequalities,  the  transition  equation,  and  ths 
switching  function  reduce  to  expressions  (49),  These  results  indicate  that 
ths  comer  point  exists  when  >  0  and  the  corner  point  1  exists 

vdisn  As  in  the  problem  of  Section  5*2,  the  two  corner  points  can- 
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(85) 
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Bodl»s  of  Cla»a  Il«  T)i«  bodies  of  this  class  are  charaotsrlssd  by 
ths  oonditlons  ($2)  trtiloh  are  valid  for 

5  *  (p  *  I  (86) 

The  geoMtry  of  the  optima  shape  can  be  obtained  froa  Eq*  (23)  and  Is 
given  by 


0  <  ?  <  1  . 


(87) 


where 


1  *  ^3^^  ^  c  ^  a  1 


(88) 


and  where 


A(11,  ♦)  ■  ± 


-  1 


♦  1  -  -  (1  -  2*1?)^ 


^  r(«,  k)  +  *(•,  k) 

2^ 


(89) 


The  syabols  F  and  E  denote  the  Incoaplete  elUptlo  Integrals  of  the  first 
and  second  kind  whose  argument  8  and  paraaeter  k  are  defined  as 
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1  -  J4  -  (1  - 


(90) 


rurth«nK>r«(  the  oonstant  f  is  relatsd  to  the  ahspe  parweeter  by  the  relation* 
ahip 


<P 


(91) 


where 


B(1)»  ♦)  •  J  [l  -  (1  -  2*Tf)^]^ 

The  drag  ooaffioleat  is  given  by 


(92) 


^  [a(0,  ♦)  T  Ad,  t)]^  [a(0,  «)  T  A(l,  ♦) 


♦  B(l,  ♦)] 


(93) 


la  the  above  expressions,  the  upper  sigis  are  valid  for  all  values  of  T) 
providing  that  1/k  *  if  g  9^,  where 


1  175) 


(94) 


otherwise,  if  cp^  <  tp  <  1/2,  the  upper  signs  are  valid  for  g  l/t/^^%  nnd  tho 
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lower  algBOi  for  t)  >  1/J^»  In  eloolngf  it  is  worth  noting  that  mom  provlous 
work  on  bodioa  of  Class  II  was  done  by  Strand  using  partly  analytical  and 
partly  nunerieal  procedures  (Ref.  4). 

Bodies  of  Class  III.  For  thisdass  of  bodies*  the  conditions  (61) 

Bust  be  satisfied  and  sure  valid  for 

1/2  <  <p  <  1  (95) 

Bence*  the  equation  of  the  optimal  shape  is  given  by 


0  <  ?  <  * 


(96) 


where  the  upper  sign  is  valid  for  £  1/Z  and  the  lower  sign,  for  t)  a  1/2. 
Finally*  the  abscissa  of  the  transition  point  and  the  shape  psuraaeter  are 
related  by  the  expression 


■  2(1  -  Ip) 

The  drag  coefficient  is  given  by 


(97) 


(98) 
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COHCLUSIOWS 

In  this  report y  the  problem  of  minimizing  the  pressure  drag  of  slender 
bodies  In  Newtonian  flow  Is  considered*  Attention  is  focused  on  the  class 
of  problems  which  Involve*  In  addition  to  the  thickness  and  the  length* 
the  enclosed  area  and  the  moment  of  Inertia  of  the  contour  for  two» 
dimensional  shapes  and  the  Mtted  area  and  the  volume  for  axlsTmmetrle  shapes. 

When  three  of  the  four  quantities  under  consideration  are  specified* 
one  obtains  a  one-parameter  family  of  optimum  shapes*  ths  parameter  being 
a  function  of  the  three  prescribed  quantities*  Qenerally  speaking*  this 
family  of  solutions  la  composed  of  three  classes  of  shapes:  (I)  a  flat 
plate  or  a  spike  followed  by  a  regular  shape*  (II)  a  regular  shape  only* 
and  (III)  a  regular  shape  followed  by  a  constant  thickness  contour  or  a 
cylinder* 

The  above  results  can  be  employed  in  order  to  understand  the  effect 
of  changing  any  one  of  the  three  given  quantities  on  the  minimum  drag  tdiile 
ths  other  two  are  kept  constant*  In  each  case*  the  drag  has  an  absolute 
miniaum*  and  the  corresponding  optimum  shapes  are  those  derived  by  Mield 
in  Refs*  1  through  3*  For  exaaqple*  if  the  thickness*  the  length,  and  the 
enclosed  area  of  a  two-dimensional  shape  are  prescribed*  the  drag  is  an  abso¬ 
lute  irfnimum  when  A/d/  ■  1/2*  which  corresponds  to  the  wedge  solution*  Kiele 
found  the  wedge  to  be  the  optimum  body  for  the  ease  where  the  thickness  and 
the  length  are  specified*  while  the  enclosed  area  and  ths  moment  of  inertia 
of  the  contour  are  free* 

In  closing*  three  remarks  are  important: 

(a)  Some  of  the  optimum  shapes  obtained  in  this  analysis  iu:e  concavei 
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oonMqtttBtly*  th«M  bodies  should  be  restudled  UBlag  the  tfewton-Busessim 
pressure  coeffioleat  law  la  secordsnoe  with  the  aethod  developed  xn  Ref*  3* 

(b)  Vftiea  the  thlokaess  ratio  beeoaes  very  large*  the  average  value  of 
the  slope  say  beooae  suoh  that  the  sleader  body  approxlnatioa  is  violated | 
eoasequeatly*  this  oase  should  be  reinvestigated  using  the  exact  tfewtoaiaa 
ejqtressioa  for  the  pressure  coefficient*  that  is*  the  sine  square  law  (Ref*  9)* 
(o)  VAmb  the  thickness  ratio  beeoaes  very  aaall*  the  frletloa  drag 
(tdiioh  was  neglected  here)  say  have  the  sane  iaportanoe  as  the  pressure  drag* 
Conaequeatly*  it  is  of  interest  to  foraulate  a  new  siiiiaal  problea  in  which 
the  total  drag  (the  am  of  the  pressure  drag  and  the  frletloa  drag)  is  ex- 
treaised  (Refs*  6  and  ?)* 
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Fig.  7.  Minumum  drag  coefficient  for  given  thickness,  length,  and  moment 
of  inertia  of  the  contour. 


Fig.  8.  Optimum  shapes  for  given  thickness,  length,  and  wetted  area. 


Fig.  9.  Shape  parameter  for  given  thickness,  length,  and  wetted  area. 
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Fig.  12.  Shape  parameter  for  given  thickness,  length,  and  volume. 


Fig.  13.  Minimum  drag  coefficient  for  given  thickness,  length,  and  volume. 


